Stability Properties of Nonhyperbolic Chaotic Attractors under Noise 
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We study local and global stability of nonhyperbolic chaotic attractors contaminated by noise. 
The former is given by the maximum distance of a noisy trajectory from the noisefree attractor, 
while the latter is provided by the minimal escape energy necessary to leave the basin of attraction, 
calculated with the Hamiltonian theory of large fluctuations. We establish the important and 
counterintuitive result that both concepts may be opposed to each other. Even when one attractor 
is globally more stable than another one, it can be locally less stable. Our results are exemplified 
with the Holmes map, for two different sets of parameter, and with a juxtaposition of the Holmes 
and the Ikeda maps. Finally, the experimental relevance of these findings is pointed out. 
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Noise plays an important role in nonlinear systems. 
Specifically, the fundamental question of the effect of 
noise on the stability of a chaotic attractor can be viewed 
under two different angles. The first aspect is to consider 
the escape from an attractor through random fluctua- 
tions. This is termed global stability. Relevant exam- 
ples range from switching in lasers , Penning traps , 
over chemical reactions kj to electronic circuits 41 . Since 
the seminal work of Kramers |5J , this problem has been 
treated for a broad range of settings [6j . For nonequilib- 
rium systems, a WKB-like extension of Kramers' equilib- 
rium theory has been devised 0,0. This so-called Hamil- 
tonian theory of large fluctuations uses an approach sim- 
ilar to path integrals, thus obtaining the most probable 
exit path (MPEP). The MPEP, with an exponentially 
favoured probability of occurrence, yields in turn the op- 
timal fluctuations and the minimal escape energy as well. 

This theory has been employed for the calculation of 
the escape from a periodic state 113 IH El Re- 
cently, it has also become possible to treat the escape 
from a nonhyperbolic chaotic attractor (NCA) [fj| , whose 
stable and unstable manifolds exhibit tangencies. It was 
demonstrated that the MPEP is uniquely determined by 
the primary homoclinic tangency (PHT) closest to the 
basin boundary. A tangency is homoclinic if both man- 
ifolds belong to the same periodic orbit and primary, if 
a perturbation is amplified under forward and backward 
iteration of the dynamics. Since, in practice, virtually all 
chaotic attractors appear to be nonhyperbolic, it can be 
considered as the general case. 

The second aspect of noise effects on NCAs is local sta- 
bility, which is a measure of the maximum distance of a 
noisy trajectory from the noisefree attractor. Here, the 
trajectory is always close to the attractor, without leav- 
ing its basin of attraction. The concept of local stability 
of a NCA against noise is of fundamental importance and 
has bearings, e.g., on noise reduction reconstruction 
of dynamical quantities [lfij| . parameter estimation [l7j . 
noise level evaluation |l8l | , and communication with chaos 



[ljfl . When applying noise bounded by a, for hyperbolic 
attractors the maximum distance scales as S max ~ a 
[20| . For NCAs, however, it was shown that there is 
a much larger b~ max as compared to the hyperbolic case 
[23, caused by attractor elonga ting deformations along 
the PHT and their images (see \TA 12.1 . as well). This 
was also confirmed experimentally [24J. 

In this Letter we contrast these two measures of sta- 
bility. While it is usually assumed that they behave in 
a similar fashion, we point here out, however, the coun- 
terintuitive effect that a nonhyberpolic chaotic attractor 
can be, in the above defined sense, globally more sta- 
ble than another one, yet locally less stable. This is 
all the more surprising as both stability properties are 
intimately related to the primary homoclinic tangency. 
This phenomenon can be understood, though, by taking 
into account that for global stability the preimages are 
most relevant, constituting the proper and uniq ue initial 
conditions for the most probable exit path [1J|. On the 
other hand, for local stability only the images govern the 
process [2l|, as their local expansion rates, given by Eq. 
(J2J below, contribute to a divergence from the attractor. 
Consequently, for local stability only linear properties of 
the system are relevant, whereas global stability can only 
be fully described by the complete set of variational equa- 
tions, which are nonlinear. 

We illustrate these findings first with the Holmes map 
[2II with two different sets of parameters. Thereafter, we 
demonstrate this phenomenon by comparing the Holmes 
and the Ikeda map 26] . Since the Hamiltonian theory 
of large fluctuations is only valid for Gaussian noise and 
the maximum distance is only well defined for bounded 
noise, we calculate for local stability also the averaged 
Gaussian distance, including higher moments. This re- 
moves any particularity of comparing different noise dis- 
tributions. The outcome of the calculation corroborates 
our main claim, too. 

As a fundamental dynamical example we consider the 
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Holmes map [2£ 



Vn+1 = a %r, 



■e. 



(1) 



with the white noise terms ^ , £ y uniformly distributed 
in the disk a: Q. + ^ < a. We choose the first set of 
parameters to be (i) a — 0.047, b — 2.4 and c = 0.155. 
That gives two attractors, symmetrical with respect to 
the origin; we focus only on one of these. When in- 
creasing 6, these attractors merge in a crisis. Our second 
set of parameters is then in the region, where only one 
large symmetric chaotic attractor exists, (ii) a = 0.01, 
b = 2.8 and c = 0.8. Both NCAs are normalized in 
a twofold way. First, the extensions in phase space 
){y max Drain 

) are demanded to be 
the same, because then the percentage of noise on each 
attractor is identical. This is a common measure of the 
relative noise intensity, in turn adjusting the local prop- 
erties. Second, the threshold of escape from the NCAs 
with bounded noise is also required to be equal. This 
guarantees the same scaling region for the maximum dis- 
tance and calibrates the global properties. For the cho- 
sen parameters, the two measures yield E « 3.1 and 



o 'escape ~ 0.13 |27j . With these two conditions met, the 
comparison is as general and unambiguous as possible. 

Let 5k be at each step of iteration the minimum dis- 
tance of the noisy trajectory from the noiseless attractor. 
The maximum distance 5 m ax of the whole trajectory is 
then defined as the maximum over all the minimum dis- 
tances: Smax — max((5fe). For the numerical computa- 

k 

tion, we partition the attractor with a grid of box edge 
length I, where I depends on the noise strength and the 
desired resolution (0.0005 < I < 0.01). We store only 
a limited number of points of the noiseless attractor per 
box of the grid (ca. 100). Each point of the noisy trajec- 
tory is then compared solely to attractor points of the box 
it falls in and the neighboring ones. If they are empty, 
the number of neighbors is increased until a point of the 
attractor has been encountered. This provides, for each 
trajectory point, the minimum distance from the attrac- 
tor 5k, and the largest of these is 5 max . With this method 
we get a much better accuracy and a larger scaling region 
than in [22I |23| , while simultaneously saving storage and 
computation time. 

The result of the calculation for the two NCAs is shown 
in Fig. ^ The scaling is limited for small noise by our 
computational resolution and for large noise by the tra- 
jectory escaping from the attractor. It is apparent from 
the graph that, for all noise intensities, set (i) (circles) 
exhibits a larger 5 rnax than set (ii) (squares), indicating 
that the attractor (ii) is locally more stable than (i). 

The two curves in the log-log scale of Fig. [I\a,re straight 
lines, interrupted by bends. Between two bends, they 
have an identical slope 1 (i.e. 5 max = P n a). The factor 
of proportionality P n varies with n, causing a different 




FIG. 1: Maximum distance S max versus noise intensity a for 
the Holmes map, with the parameters (i) a — 0.047, b — 2.4 
and c = 0.155 (circles) and (ii) a = 0.01, b = 2.8 and c = 0.8 
(squares). The regions of linear growth are fitted by straight 
lines and marked with the number n of the corresponding im- 
age of the PHT. For each noise strength 5 x 10 9 iterations have 
been used. The inset shows the average (5 q ) 1 ^ q , q = 1, 10, 
with Gaussian noise for set (i) (circles) and set (ii) (squares). 
10 9 iterations for each noise level have been averaged over. 



offset. The 5k achieve their maxima at the PHT and 
images thereof (see Fig. 5 of [2lJ for a very instructive 
illustration) . At the n — th image, a perturbation at the 
PHT grows like ^ 



+ 1, (2) 



where /(x) is the dynamics of the system at x, Df(x.) 
the Jacobian, and e™_f ^ [/■? (x)] the most expanding unit 
vector at / J '(x) under the application of D/ n ~ ! (x). This 
factor sums up all the maximal stretching factors of the 
PHT and its images up to the n — th one. Typically, 
Pn > P-m for n > m, implying that the distance from the 
attractor increases with the number of images. However, 
because higher images are folded back on the attractor, 
other parts of the attractor instead of an iteration of the 
PHT come closer to the noisy trajectory as the noise level 
is incremented. This results in a saturation of the max- 
imum distance, which produces a bend. In turn, when 
the noise strength is further augmented, 5 max switches 
to the next lower image of the PHT, again initiating a 
regime of linear growth, and so on. 

In Fig. ^ the sequences 4^3-^2 for set (i) and 
2 — > 1 for set (ii) can be seen. For set (ii), the offsets 
from the numerics of Fig. ^ for n = 1,2 are 8.3,9.5 
(solid lines), while Eq. J2J yields P n — 8.2,9.5, a very 
good agreement. Set (i) for n = 2,3,4 gives 15,17,34 
from Fig. 2] (solid lines), whereas Eq. Q results in 
P n = 13.4, 16, 34, also a reasonably good agreement. The 
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values for lower images of the PHT (i. e. higher noise) fit 
slightly worse. However, the matching can be improved 
by using the full dynamics instead of the linearized Eq. 
PJl. since nonlinear effects play an increasing role for 
larger noise levels. By doing this, one gets 14.5, 16.5, 34, 
again in good accordance. 

To provide a better basis for the comparison with 
global stability, we calculate the averaged moments of 
the distance = (■£ J2k=i($k) q ) 1/q using Gaussian 

white noise, with = and = cr 2 8ij. This is 

shown in Fig. ^ inset, for q = 1, 10. The correspond- 
ing moments for set (i) are for all q above the ones of 
set (ii), more distinctive for higher q. The same applies 
for (5 q ) 1 / q with bounded noise (not shown). Here, in 
the limit q — > oo the maximum distance is recovered 

\ u / 7 u max- 

Global stability is evaluated with the Hamiltonian the- 
ory of large fluctuations, solving a variational equation 
for the MPEP [H HE!, which provides the action 
S = \ Sn=i ^n, with A n the optimal fluctuations. 
The mean first exit time is then given by (r) ~ exp [-^-1 . 
The MPEP starts at the preimages of the PHT, leaves 
the attractor close to the PHT and moves along their im- 
ages towards the saddle point on the basin boundary ^3 • 
Employing this scheme, one obtains for set (i) S « 0.015 
and for set (ii) S « 0.01, meaning now that set (i) is glob- 
ally more stable than set (ii). We stress that this leads, 
e.g. for a noise value of a 1 = 0.001, to an amplification 
of (t) by a factor of exp [ 0, °2 5 ] rj 148, it is therefore no 
small effect. 

These opposing stability properties establish our main 
result. The Holmes map (as a typical example of a NCA) 
is with set (i) of parameters locally less stable than with 
set (ii), i.e., the maximum distance 8 ma x is larger, but 
globally more stable, i.e., the escape energy and conse- 
quently the mean first exit time are larger. 

Next we demonstrate that this phenomenon can be 
much more pronounced when comparing two NCAs orig- 
inating from different dynamical s yste ms. For that pur- 
pose, we introduce the Ikeda map [26j 



Zn+i = a + bz n exp 



1 



(3) 



where z n = x n + iy n . We fix the parameters at a = 
0.9, b — 0.9, k — 0.3 and rj = 6.0, which results in a 
NCA. We compare this NCA with the one obtained for 
the Holmes map with the parameter set (iii) a = 0.01, 
b = 2.78 and c = 1.56. Again both attractors are normal- 
ized in the two ways explained above, with E ps 2.1 and 
Cescape ~ 0.1 . The maximum distance is depicted in 
Fig. The features are more striking than in the pre- 
vious example, 8 max differs, for instance, for a = 10~ 4 , 
by one order of magnitude. Furthermore, the scenario of 
jumping from one image of the PHT to the next one hap- 
pens for the Ikeda map more frequently. For the lowest 
noise level considered, the maximum distance occurs at 




FIG. 2: Maximum distance 8 max versus noise level a for the 
Ikeda and the Holmes maps, with the parameters (Ikeda) a = 
0.9, b = 0.9, k = 0.3, and r\ — 6.0 (circles) and set (iii) a = 
0.01, b — 2.78 and c = 1.56 (squares). For each noise strength 
5 x 10 9 iterations have been used. The inset shows the average 
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g = l, 10, with Gaussian noise for Ikeda (circles) and 



set (iii) (squares) . 10 9 iterations for each noise level have been 
averaged over. 



the 6-th image of the PHT. 

For set (iii) of the Holmes map, the numerical off- 
sets of Fig. |21 (solid lines) come about as 8.2,11 for 
the images n — 1,2 of the PHT. Equation re- 
sults in P n = 8, 11.25, agreeing extremely well. The 
Ikeda map gives for n = 2,3,4,5,6 the numerical val- 
ues 10, 33, 53, 95, 225, respectively, while evaluated with 
Eq. @ produces P n = 8,22,50,96,247. Taken into ac- 
count that on the one hand for large images of the PHT 
the maximal distance is numerically hard to observe, as 
several subsequent optimal fluctuations are needed to 
achieve it, and on the other hand for low images the 
noise is already so large as to cause nonlinear effects, the 
correspondence is tolerably good. 

In the inset of Fig. [3 the averaged distances with 
Gaussian noise {5 q ) 1 ^ q ,q = 1, 10 are displayed. For small 
q the Holmes map is here above the Ikeda map. This 
is rooted in the fact that the unstable manifold of the 
Ikeda map is more curved at the PHT and their images. 
Hence, the average exhibits less of the maximal possible 
expansion. However, for larger q (g = 10 in the graph), 
the average is above for all noise levels. Again, the same 
holds for (S q ) 1 ^ q and bounded noise (not shown). 

Global stability analysis, as before, entails for the Ikeda 
map S fss 0.025 and for the parameter set (iii) of the 
Holmes map S ~ 0.007. Consequently, for the selected 
parameters, the Ikeda map is globally much more stable 
than the Holmes map, while it is locally much less sta- 
ble, which is caused by the higher images of the PHT 
having larger expansion factors [Eq. @] and at the 
same time weaker folding back to the attractor. Both 
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effects are most pronounced in the reievant low noise 
limit. The Ikeda map is globally more stable by a factor 
of exp r 0,0 j- 8 ] . The amplification becomes huge for small 
noise (e.g. 6.5 x 10 7 for a 2 = 0.001) and is easily measur- 
able. This establishes that the phenomenon of opposite 
stability properties, when comparing NCAs originating 
from different dynamical models, can be observed in an 
even more striking manner. 

We have confirmed this counterintuitive phenomenon 
also when comparing the Henon map with both, the 
Ikeda and the Holmes maps, corroborating our findings, 
which we claim to be a general feature of NCAs. 

In the present work, we were not concerned with the 
overall scaling of S max , only with the fact that one curve 
lies above another, thus implying being locally less sta- 
ble. In [2^|. however, it was claimed that the scaling is 
Smax ~ cr 7 , were 7 = 1/D±, with D\ the information 
dimension of the attractor. The agreement between this 
value and our, very accurate, numerics, is not too good, 
though [23. This discrepancy is caused by the fact that 
in the derivation of the scaling in not P n from Eq. 
@ was used, but the positive Lyapunov exponent, which 
usually has a smaller value. Thus, in general 1/D\ can 
be regarded only as a lower bound for 7. The question 
of an exact scaling exponent will be treated in |30j | . 

As our findings can have an huge effect on the max- 
imum distance and the average escape time, they have 
also relevance for experiments, since one cannot simply 
and straightforwardly conclude the behavior of one of the 
stability types by measuring the other. 
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